cd D:\北工大\B503\合作论文\注册制和金融效率\code\stata
* step1
clear
import excel "D:\北工大\B503\合作论文\注册制和金融效率\数据\通过检验\psm_data_2019-8.xlsx", sheet("Sheet1") firstrow

//将省份转变为可回归的
encode 省份 ,gen(prov)
//行业也是提取第一个字母作为行业分类
tab 行业级别
gen industrial=substr(行业级别,1,1)
tab industrial
encode industrial ,gen(industry)
//reg y i.treat##i.D x1 x2 x3 x4 x5 x6  
//创建交乘项
gen did=treat*D

keep if inlist(行业级别a, 1)

xtset id Time 
//进行did回归。可以在x5 后面加入其他控制变量，以及加上i.id或者别的进行控制
// xtreg y did x1 x2 x3 x4 x5 x6 ,fe be re

xtreg y did x1 x2 x3 x4 x5 x6, re
est store r1
esttab r1  using 行业级别0异质性分析.doc, b(%12.3f) se(%12.3f) nogap compress replace s(N r2 r2_a r2_w r2_b r2_o ) star(* 0.1 ** 0.05 *** 0.01)





//
// local chi2value = 14.81  // 假设的 chi2 值  
// local pvalue = 0.0386   // 假设的 Prob > chi2 值  

//内生性检验
// // Hausman
// xtreg y did x1 x2 x3 x4 x5 x6 x7 , re
// estimates store ols
// // x1是内生变量，其工具变量是Z1和Z2（以及可能的其他未列出的工具变量） 
// // x2至x3（以及可能的其他未列出的变量）作为控制变量
// ivregress 2sls y (x5 = x6 ) x1 x2 x3 x4  x7 //两阶段最小二乘法
// // ivregress 2sls y (x1 = x5 x6)  x2 x3 x4  x7 //两阶段最小二乘法
//
// estimates store iv //存储回归结果
// outreg2 using "my_regression_results.doc", word  dec(3)  ctitle(iv) replace
// hausman iv ols, constant sigmamore //Hausman检验
// outreg2 using "xtreg", excel word dec(3) adds(chi2, r(chi2), Prob > chi2, `r(p)') adec(3) ctitle(hausman) append
// estat endogenous   // DWH检验
// estat firststage // 弱相关性检验


